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Abstract. We prove a large deviation principle of Freidlin-Wentzell's type for the multi- 
p I ' valued stochastic differential equations with monotone drifts, which in particular contains 

—1 a class of SDEs with reflection in a convex domain. 
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1. Introduction 

Consider the following multivalued stochastic differential equation (MSDE in short): 



f dX(t) G b(X(t))dt + a(X(t))dW(t) - A(X(t))dt, 
O i I X(0) = x G D(A), 1 > 

where A is a multivalued maximal monotone operator, which will be described below, 
W(t) = {W k (t),t ^ 0, k G N} is a sequence of independent standard Brownian motions 
on a filtered probability space (fi, J 7 , P; (J*t)^ Q ), 6 : M m -»■ M m and a : M m ->> M m x / 2 are 
two continuous functions, / 2 stands for the Hilbert space of square summable sequences 
of real numbers. 

^ ! This type of MSDE was first studied by Cepa in [5j E]. He proved that if b and a are 

Lipschitz continuous, then there exists a unique pair of processes (X(t),K(t)) such that 

X{t) = x+ [ b{X{s))ds+ [ a(X(s))dW(s)-K(t), 
Jo Jo 

where K(t) is a process of finite variation (see Definition 12.31 below for more details). 
Recently, Zhang p3] extended Cepa's result to the infinite dimensional case, and relaxed 
the Lipschitz assumption on b to the monotone case. It should be noted that when 
A is the sub differential of the indicator function of a convex subset of R m , the above 
MSDE is the same as the usual SDE with reflecting boundary in a convex domain (cf. 
[H,[n]). Moreover, since the sub differential of any lower semicontinuous convex function is 
a maximal monotone operator, Cepa's result can also be used to deal with the SDE with 
discontinuous coefficients. It is well known that there are many literatures to investigate 
the SDEs with reflecting boundary since the solutions of a large class of PDEs with 
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Neumann boundary and mixed boundary conditions can be represented by the solution 
of such SDEs (cf. pP). 

We now consider the following small perturbation of Eq.flTJ): 

f dX £ (t) G b( X e (t))d t + y/Za(X e (t))dW(t) - A(X e (t))dt, 

\ X e (0) = x G D(A), e G (0,1]. { ] 

The solution of this equation is denoted by (X e (t, x), K e (t, x)). We want to establish the 
large deviation principle of the law of X e (t,x) in the space § := C([0,T] x D(A); D(A)), 
namely, the asymptotic estimates of probabilities P(X e G T), where T G B(E>). 

In [I], Anderson and Orey considered the same small random perturbation for the 
dynamical system with reflecting boundary in smooth domain, and obtained the Freidlin- 
Wentzell's large deviation estimates in C([0, T]; D(A)). They assumed that the coefficients 
are bounded and Lipschitz continuous, and the diffusion coefficient is non-degenerate. 
Using the contraction principle, Cepa [5] only considered the large deviation principle 
of one dimensional case based on an explicit construction of the solution (cf.[13]). The 
mult i- dimensional case is still open. Compared with the usual SDE, i.e., A = 0, most 
of the difficulties come from the presence of the process of finite variation, K(t). One 
only knows that 1 1— > K(t) is continuous, and could not prove any further regularity such 
as Holder continuity. Therefore, the classical method of time discretization is almost 
inapplicable (cf. [7]). 

Our method is based on the recently well developed weak convergence approach due to 
Dupuis and Ellis [8] (see also [2j [3]). This method has been proved to be very effective 
for various systems (cf. [TTJ [T5J HJ EH E3HJ [121 etc.]). In the situation considered in the 
present paper, however, since we cannot prove the following uniform estimate as in [TT] : 
for any p ^ 2 and s,t G [0, T], x,y G D(A) 

sup E\X € {t,x) -X e {s,y)\ 2p ^ C{\t-s\ p + \x-y\ 2p ), 

eS(0,l) 

we cannot obtain the tightness of the laws of X e (t, x) in 8. Some technical difficulties for 
verifying the conditions (LD)i and (LD)2 below need to be overcome. 

In Section 2, we recall some well known facts about the MSDE and a criterion for 
Laplace principle. In Section 3, we present our main result and give a detailed proof. 
Throughout the paper, C with or without indexes will denote different constants (de- 
pending on the indexes) whose values are not important. 



2. Preliminaries 

We first give some notions and notations about multivalued operators. Let 2 Rm be the 
set of all subsets of MJ 71 . A map A : W 71 — > 2 Rm is called a multivalued operator. Given 
such a multivalued operator A, define: 

D(A) := {xeR m : A(x) ^ 0}, 
lm(A) := U xeD{A) A(x), 
Gt(A) := {(x, y) G M. 2m : x G R m , y G A{x)}. 

We recall the following definitions. 

Definition 2.1. (1) A multivalued operator A is called monotone if 

(yi - 2/2, xi - x 2 ) Rm ^ 0, V(xi, yi), (x 2 , y 2 ) €. Gr(A). 



(2) A monotone operator A is called maximal monotone if for each (x,y) G Gr(A), 

(y - y',x -x') Rm ^ o, \/(x',y') eGx(A). 

Examples 2.2. Suppose that O is a closed convex subset ofW 1 , and Iq is the indicator 
function of O, i.e, 

r o, if x g o, 
w - 1 +00; if x ^ a 

The subdifferential of Iq is given by 

dl (x) = {yeW n :(y,x-z) Rm ^0,\/zeO} 

0, if x i O, 
{0}, if x G Int(O), 
A x , if x G dO, 

where Int(O) is the interior of O and A x is the exterior normal cone at x. One can check 
that OIq is a multivalued maximal monotone operator in the sense of Definition \2. 1[ 

We now give the precise definition of the solution to Eq.flT]). 

Definition 2.3. A pair of continuous and (J~t)- adapted processes (X,K) is called a so- 
lution of Eq. (1) if 

(i) X(0) = x, and for all t ^ 0, X(t) G D(A) a.s.; 

(ii) K(0) = a.s. and K is of finite variation; 

(ui) dX(t) = b(X(t))dt + a(X(t))dW(t) - dK(t), < t < oo, a.s.; 
(iv) for any continuous and (J~t)- adapted processes (at, (3) with 

(a(t), f3(t)) G Gt(A), Vt G [0, +oo), 

the measure 

(X(t)-a(t),dK(t)-P(t)db)^0 a.s.. 

We now recall an abstract criterion for Laplace principle, which is equivalent to the 
large deviation principle (cf. 0,m[T5]). It is we U known that there exists a Hilbert space 
so that I 2 C U is Hilbert-Schmidt with embedding operator J and {W (t),k G N} is a 
Brownian motion with values in U, whose covariance operator is given by Q = Jo J*. For 
example, one can take U as the completion of I 2 with respect to the norm generated by 
the scalar product 



For a Polish space B, we denote by B(E) its Borel a- field, and by Cy(B) the continuous 
function space from [0, T] to B, which is endowed with the uniform distance so that Qr(B) 
is still a Polish space. Define 



T 



h(s)ds : heL 2 (0,T-l 2 ) } (3) 



with the norm 

/ rT ^ V2 

\\hy T ■= ( / \\h(s)\\ 2 2 ds 



where the dot denotes the generalized derivative. Let \x be the law of the Brownian motion 
W in C T (U). Then 

(4,c T (u),/i) 
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forms an abstract Wiener space. 
For T,N > 0, set 

V N := {h e f T : \\h\\%^N} 

and 

^ T j h : [0, T) — y I 2 is a continuous and (J-i)-adapted ^ 
A 1 process, and for almost all u, h(-,u) G V N 

We equip T>^ with the weak convergence topology in ij,. Then 

"Div is metrizable as a compact Polish space. (5) 

Let § be a Polish space. A function I : § — y [0, oo] is given. 

Definition 2.4. TTie function I is called a rate function if for every a < oo, the set 
{/ G § : J(/) ^ a} is compact in S. 

Let {Z e : Ct(U) — )■ S, e G (0, 1)} be a family of measurable mappings. Assume that 
there is a measurable map Zq : l\ i— )■ § such that 

(LD)i For any iV > 0, if a family {/i e , e G (0,1)} C (as random variables in T>n) 
converges in distribution to h G «4^-, then for some subsequence e^, Z efc ' 



converges in distribution to Z (h) in S. 

(LD) 2 For any iV > 0, if {h n ,n G N} C T>n weakly converges to h G then for some 
subsequence h nh , Z (h nk ) converges to Z (h) in §. 

For each / G S, define 

1(f) '=7; inf NI/>> (6) 

where inf = 00 by convention. Then under (LD)2, /(/) is a rate function. 
We recall the following result due to [3] (see also [T71 Theorem 4.4]). 

Theorem 2.5. Under (LD)! and (LD) 2 , {Z e ,e G (0, 1)} satisfies the Laplace principle 
with the rate function 1(f) given by (0|). More precisely, for each real bounded continuous 
function g on E>: 



lime log E M exp 



g(z e 



-mfRf) + /(/)}. (7) 

JtSs 



In particular, the family of {Z e ,e G (0,1)} satisfies the large deviation principle in 
(§,£?(§)) u>rf/i i/ie raie function 1(f). More precisely, let v e be the law of Z e in (§,£?(§)), 
t/ien /or any i? G 

— inf /(/) ^ liminf elogvAB) ^ limsupelog^ e (i?) ^ — inf /(/), 

feB° e->0 e ^ / e _g 

where the closure and the interior are taken in S, and 1(f) is defined by (f/^j. 

3. Main Result and Proof 

We assume that 

(HI) A is a maximal monotone operator with non-empty interior, i.e., lnt(D(A)) ^ 0; 
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(H2) a and b are continuous functions and satisfy that for some C a , Cb > and all 

x, y G M m 

- 0-(y)|| L2( ,2. Rm) < C a \x-y\, 
(x - y,b(x) - b(y)) Rm ^ C b \x-y\ 2 , 
where L 2 (/ 2 ;M m ) denotes the Hilbert-Schmidt space and | • | denotes the norm in 



i m , and for some C' b > and n G N 



i^)Ka(i + |xr). 

It is well known that under (HI) and (H2), there exists a unique solution (X e , K e ) to 
Eq.(T5]) in the sense of Definition 12.31 (cf. [H]). Our main result is stated as follows: 

Theorem 3.1. Assume that (HI) and (H2) hold. Then the family of {X e (t,x), e G 
(0, 1)} satisfies the large deviation principle in § := C([0, T] x D(A); D(A)) with the rate 
function given by 

2 inf JWlL ( 8 ) 

where X h (t,x) solves the following equation: 

dX h {t) G b{X h {t))dt + a{X h {t))h{t)dt - A(X h (t))dt, X h {0) = x. 



For proving this result, by Theorem 12.51 the main task is to verify (LD)i and (LD)2 
with 

§:=C([0,T] xD(Aj;D(Aj), Z e = X e , Z (h) = X h . 
This will be done in Lemmas 13.71 and 13.81 below. 

Let h e G A% converge almost surely to h G Ajj as random variables in an d 
(X £ ' h % K e ' e ) solve the following control equation: 



X e ^(t) = x+ b(X e ' he (s))ds + / (x(X eA (s))/i e (s)ds 
Jo Jo 

f a{X e ' h ^{s))dW{s) - K e ' h ^{t), (9) 
Jo 

which can be solved by Girsanov's theorem, and (X h , K ) solve the following deterministic 
equation: 



X h (t) = x + [ b(X h (s))ds+ [ a(X h (s))h(s)ds - K h (t). 
Jo Jo 



(10) 



Let \K\ s t denote the total variation of K on [s,t\. We recall the following result due to 
Cepa [6] (see also [MJ Propositions 3.3 and 3.4]). 

Proposition 3.2. Under (Hl) ; there exist a G lR m ,7 > 0,/i ^ such that for any pair 
of (X, K) with the property (iv) of Definition \2.3\ and all ^ s < t ^ T 

(X(r) - a, dK(r)) Rm ^ ^y\K\ s t - // J \X(r) - a\dr - lt i{t - s). (11) 

Moreover, for any pairs of (X,K) and (X,K) with the property (iv) of Definition \2.?A 

(X(t) - X(t), dK(t) - dK(t)) Rm > 0. (12) 

Using this property, we first prove the following uniform estimates. 
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Lemma 3.3. For any p ^ 1, there exists C P) t,n > such that for any e G (0, 1) and 
x,ye D(A) 



E( sup \X*> h %x)-X*> h %y)\ 2p ) ^C p ^ N \x-y\ 2p . 
,te[o,T\ 



(13) 



Proof. Set 
and 



Z e (t) :=X*> h *{t,x)-X^{t,y) 



A(s) :=a(X^(s,x))-a(X^(s,y)). 
By Ito's formula, (H2) and ffl2l) . we have for any p ^ 1 

|Z £ (t)|^ = \Z e (0)\ 2p + 2p f \Z e {s)\^- 2 {Z e {s)XX^{s,x)) -b(X^(s,y))) Rm ds 

Jo 



+2p / \Z e (s)\ 2p - 2 (Z e (s),A(s)h e (s)) Rrn ds 
Jo 

+2pV~e [ t \Z e (s)\ 2p - 2 (Z e (s),A(s)dW(s)) K 



-2p C \Z € (s)\ 2p - 2 (Z € (s),dK^(s,x) -dK^^y)}, 
Jo 



-p / |Z e ( S )|^(||A( S )r + 2( P -l)(Z e ( S ),A( s )A*( S )Z e ( s )) Km /|Z £ ( S )|^)d S 
Jo 

^ \Z e (0)\ 2p + C f \Z e {s)\ 2p ds + 2p f \Z e (s)\ 2p ■ \\h e (s)\\ P ds 
Jo Jo 

+2pV~e [ t \Z e (s)\ 2p - 2 (Z e (s),A(s)dW(s)) Rm . 
Jo 



Set 



g(t):=E\ sup \Z e {s)\ 2p 

K se[o,t] 

By BDG's inequality and Young's inequality, we have for any 5 > 



E 



sup / \Z e (s)\ 2p ~ 2 (Z e (s),A(s)dW(s)) R 
t'e[o,t] Jo 

^CE fj\z e (s)\^A(srZ e (s)\\^ 

< CE ( sup \Z e {s)\ 2p [ \Z e {s)\ 2p ds 
\se[o,t] Jo 

^5-g(t) + C s f E\Z e (s)\ 2p ds. 
Jo 



1/2 



(14) 



Similarly, we have 



E 



\Z e (s)\ 2p ■ \\k(s)\\ l2 ds 



^ NE \Z e (s)\ 4p ds 



1/2 



^ 5-g(t) + C s ,N f nz t (s)\ 2v &s. 
Jo 



(15) 



Letting 5 = 1/4 in (lT4"j) and f|T5|) and combining the above calculations, we get 



t rt 



9(t)<g(0) + lg(t)+C I E\Z e ( S )\ 2 ?ds ^ 2g(0) + 2C I g{s)ds, 
1 Jo Jo 

which gives the desired estimate by Gronwall's inequality. □ 



Lemma 3.4. For any p ^ 1 and x G D(A), there exists C Pj t,n,x > such that for any 
e e [0, 1) 

E ( sup \X^(t,x)\ 2p ) +E\K e ' h *(-,x)\ T ^ C PjT)N , x . (16) 
\te[o,T] J 

Proof. First of all, as in the proof of Lemma 13.31 we can prove that 

E[ sup \X^(t,x)\ 2 A ^C P>T>N , X . (17) 
\te[o,T] J 

Let a G lnt(D(A)) be as in Proposition 13.21 By Ito's formula, ( ITT1) and (H2), we have 
V eA (0-«| 2 = l\x-a\ 2 + I' (X^(s)-a,b(X^(s))} Rm ds 



2' w 2 

+ / (X^(s) - a,a(X^(s))h e (u)) Rm ds 



Jo 

t 







/ (X e ^{s) -a,a(X^(s))d^(s))j 
Jo 

- f (X^(s)-a,dK^(s)) Rm 
Jo 



+ 5 / \\a(X^(s))\\ 2 ^ds 



2 



Jo 



L 2 {l 2 ;M. m P 



^ -\x-a\ 2 + C b \X e ' h *(s) -a\ 2 ds 
2 Jo 



+ / (X^(s) -a,6(a)) Rm ds 
Jo 

+iV Qf | ( j(X e ' /le (s))*(X e ' /! ' e (s) - a)| 2 d^ 

/ (X e ^(s) - a, a(X e ^(s))dJy(u)) Rm + /z 7 t 

- 7 |^f + A* [ \X^(s)-a\ds 
Jo 

[\\X^( S )\ + a(0)) 2 ds. 
Jo 



+ Cle " 



The desired estimate now follows by ( FlTl) . □ 
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Define 

-hi 



w e (t,x):= / a(X h (s,x))(h e (s) -h(s))ds. (It 
Jo 

The following lemma is easy by Ascoli-Arzela's lemma. 



Lemma 3.5. w e (-,x) converges a.s. to zero in C([0,T], D(A)) . 

We now prove the following key lemma. 
Lemma 3.6. X e,h6 defined by (OJ) converges in probability to X h defined by ^(%) in S. 
Proof. Set v e (t) := v e (t,x) := X e ' h *(t,x) - X h (t,x). Then 

v e (t) = K e > h '(t)-K h (t)+ [\b(X e > h *(s))-b(X h (s)))ds 



+ / {a{X e > h *( S ))h e {s)-a(X h (s))h(s))ds 
Jo 

+V~e [ o-(X e > h *(s))dW(s). 
Jo 

By Ito's formula, we have 

\v e (t)\ 2 = 2 / (v f (s),dK^(s)-dK h (s)) Rtn 
Jo 

+2 / (v € (s),b(X^(s))-b(X h (s))) Rm ds 
Jo 

+2 / (v £ (s),(a(X*> h <(s)) - a(X h (s)))h e (s)) Rm ds 
Jo 

+2 f (v e (s),a(X h (s))(h £ (s) - h(s))) Rm ds 
Jo 

+2V~e [ (v e (s),a(X^(s))dW(s)) Rm 
Jo 

\\a(X^(s))\\l 2{mm) ds 

= : ij(t) + / 2 e (t) + / 3 e (t) + ii(t) + im + mt). 
m < o 



It is clear that by (112 
and 



r 2 (t) ^ 2C b [ 

Jo 

By BDG's inequality and (H2) we also have 



t 

2, 



\v e (s)\ 2 ds. 



E sup |J 5 e (t)| + E sup ) <CVi. 

,te[o,T] / Vte[o,T] 



As estimating f|T5|) we have 



E ( sup m(s)\ J < ^E ( sup |t» e (s)| 2 J +C / E|t» e (s)| 2 ds. 
Vse[o,t] y 2 Vse[o,t] y Jo 



Set 



Then we have 



g(t) := E [ sup \v e (s) 
,«e[o,t] 



which implies that 



9(t) < lg(t) + C^ + E\ sup |J|(s)| ) + C fnv^lMs, 
1 \se[o,t] J Jo 

g{t) «C (Vi +2E ( sup |/|(s)|]+C / </(s)ds. 

\s6[0,T] y Jo 

By Gronwall's inequality we get 

E( sup \v e (s)\ 2 ) ^C^+CE\ sup (19) 
\-e[o,n y \se[o,T] J 

We now deal with the hard term J|. By Ito's formula again, we have 

1 



2 r 4 (t) = (v e (t),w e (t))^- J o w e (s)d(K^(s)-K h (s)) 
- f w e (s)(b(X^(s)) -b(X h (s)))ds 







- / w e (s)(a(X e '^(s))h e (s) - a(X h (s))h(s))ds 



o 



-yfe J w e (s)a(X e ^{s))dW{s) 

=: F 41 (t) + Ih(t) + lUt) + IIS) + IIS)- 
For any 5 > and R > 0, we have 

P ( sup |J^(t)| ^ 5 ) = P ( sup ^ 5; sup Mi)) < P 

\te[o,T] J \te[o,T] t<=[0,T] 

+P ( sup \ItS)\ > 8; sup |u e (*)| ^ P 
Vte[o,T] te[o,T] 

^ P I sup |w e (t)| ^ 5/P J + W sup |u e (*)| ^ P 
ye[o,T] y \te[o,T] 

By Lemma 13.51 and (1161) we know 

limP f sup ^ cH =0. 

\<e[o,T] y 

Noting that 

sup \It 2 (t)\ ^ sup | We ( a )| • {\K^{-)\%+ \K h {-)\%). 
te[o,T] se[o,T] 

as above, we also have 

sup (/^(t)! — > in probability. 

te[o,r] 
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Similarly, we have 



sup |/43(i)|+ sup I /| 4 (2) I — )■ in probability. 

*e[o,T] *e[o,T] 



Moreover, by BDG's inequality we have 



V^E [ sup |/ 4 e 5 (2)| ) < C^~e. 
te[o,T] 



Combining the above calculations, we get 



sup ^(i) I — > in probability. 

te[o,T] 



It is easy to see by (1161) that 



Hence 



e->0 

Substituting this into f|T9|) we obtain 



sup E sup ^(t) | 2 ] < +00. 
ee(o,i) Vte[o,T] 



limE ( sup |/ 4 e (2)| ] = 0. 
te[o,T] 



limE sup \v e {s)\ 2 ) = 0. 
\se[o,T] 



Thus, we have proven that for all x G D(A) 



sup \v e (t, x)\ 2 — > 0, in probability. 

te[o,T] 

We now strengthen it by Lemma 13.31 to 

£ n>e := sup |^e(2, x)\ 2 — > 0, in probability. 

te[0,T\,xeD(A),\x\^n 

Set 

:= D(A) n{i£K m : \x\ < n} n <5Z m , 

where 5 > and <5Z m denotes the grid in M. m with edge length 5. It is clear that there 
are only finite many points in 3 s n . Hence 

£* e := SU P \ v e(t, x)\ 2 — y 0, in probability. 

te[o,T],xev>i 

For any a; G M m , let denote the left-lower corner point in 5Z m so that 

\% — xg\ ^ 8- 

Noting that 

£ n)£ ^2£* e + 2 sup sup \v e (t,x) - v € (t,x s )\ 2 , 

we have for any > and some a > 

P(^,e>4/3) ^ P(&, e + _sup sup |v e (t,a;) -^(2,x <5 )| 2 ^ 2/3) 

xeD(A),|x|s:n*6[0,r] 

< P(e^>/3) + Pl _sup sup \v e (t,x) -v e (t,x s )\ 2 > (3 

\xeD{A),\x\^nt£[0,T] 
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< P(&,e> _SUP SUp \v e {t, X ) -V e {t,X S )\ 2 ]/P 

\x<ED(A),\x\^nte[0,T] J 

where the last step is due to Lemma 13.31 and Kolmogorov's criterion. First letting 5 be 
small enough, then e go to zero, we then obtain 

limP(£„, £ >4/3) = 0. 

which yields the desired result. □ 
Lemma 3.7. (LD)i holds. 

Proof. Let h e be a sequence in Aj? converge to h in distribution. Since is compact 
and the law of W is tight, {h e , W} is tight in T>n x Ct(U) by the definition of tightness. 
Without loss of generality, we assume the law of {h e , W} weakly converges to /x. Then 
the law of h is just /i(-,Ct(U)). Indeed, for any bounded continuous function g on IV, 
we have 

E(g(h)) = lim E(g(h e )) = [ g(h)fi(dh,C T (V)). 

J V N 

By Skorohod's representation theorem, there are (O, P) and {h e , W e } and {h, W} such 
that 

(1) (h t ,W £ ) a.s. converges to (h,W); 

(2) (h e , W £ ) has the same law as (h e , W)\ 

(3) The law of {h, W} is /i, and the law of h is the same as h. 
Using Lemma [3. 6 [ we get 

$(-r= / h e ds + W e ) -> X~ h , in probability, 
v e Jo 

where $ is the strong solution functional (cf. [9]). From this, we derive 

$(— t= / h e ds + W) -> X h , in distribution. 
v e Jo 

Thus, (LD)i holds. □ 

Similar to the proof of Lemma I3.6[ one can easily verify that 
Lemma 3.8. (LD) 2 holds. 

Thus, by Lemmas 13.71 13.81 and Theorem 12.51 we have proved Theorem 13.11 
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